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This  paper  characterizes  a  class  of  constrained  transportation 
problems  (linear  programming  problems  which  are  composed  of  a 
transportation  problem  with  additional  constraints)  which  can  be  trans¬ 
formed  into  equivalent  transportation  problems.  An  efficient  computational 
procedure  for  finding  an  equivalent  transportation  problem  is  presented 
which  introduces  at  most  one  origin  and  one  destination  for  each  extra 
constraint.  Our  results  extend  the  rrocedures  developed  by  Wagner  and 
Manne  for  problems  in  which  the  extra  constraints  consist  of  bounding 
certain  partial  sume  of  variables. 


DD  ?0~.1473  (PAGE  1) 


s  'k  oi  oi  -Acr-oei : 


Unclassii  <ed 


DD  '~..1473  =«*» 


1.  :xt?.opjctio* 


Linear  progr&rcing  literature  abounds  with  ingenious  procedures  for 
reformulating  linear  problems  wnicn  seer  unrelated  to  tne  transportation 
orjt'ter  into  ar,  equivalent  transportation  problem  [  2,3,4,5,10,11,12]. 

Sore  of  these  linear  probiers  are  constrained  transportation  profci ers 
'7,8]  (i.e.,  a  linear  progranring  probler  consisting  of  a  transportation 
probier  with  extra  constraints).  In  this  paper  we  substantially  enlarge 
the  class  of  constrained  transportation  probier  formulations  which  nave 
Peer,  snowr.  to  be  transportable  into  an  equivalent  ordinary  t.  cnsportation 
probler  by  extending  tne  procedures  developed*  by  wagner  [10]  and  Hanne 
[5,  p.  332].  There  are  decided  computational  advantages  to  transforming 
a  seeringiy  mo re  general  linear  programing  probler  into  a  transportation 
probler  when  this  is  possible.  For  instance,  the  computational  results  of 
]1,6,9]  disclose  tnat  transportation  problems  can  be  solved  at  least  100 
ti-es  faster  oy  spec  al  purpose  codes  tnan  by  general  purpose  linear  pro- 
granriog  -etnocs.  Trus,  trve  difference  oetween  soiling  a  constrained  trans- 
bo-taticr  probler  witn  a  general  setnod  and  solving  its  counterpart  trans¬ 
portation  probler  witn  a  special  purpose  code  should  oe  very  substantial, 
certainly  not  less  tnan  the  factor  of  100  to  1. 

Tnere  are,  of  course,  other  advantages  to  transforming  a  constrained 
transportation  probler  into  a  transportation  probler.  By  using  the  special¬ 
ized  algoritrrs  [1,6,9],  no  round-off  error  is  introduced.  Also,  if 
tne  total  supply  of  the  reformulated  problem  is  not  equal  to  tne  total  de- 
~anc,  tnen  cne  can  conclude  without  further  effort  tnat  tne  original  prob¬ 
ler  lacrs  a  feasible  solution. 

«agner  [10]  ana  *<anne  [5,  p.  332]  showed  tnat  if  an  extra  constraint 
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to  a  transportation  problem  bounds  a  partial  sum  of  variables  in  a  single 
node  constraint,  then  the  problem  can  be  reformulated  as  a  larger  trans¬ 
portation  or  network  problem  with  addition!  nodes.  This  paper  extends 
the  tecnniques  proposed  by  Vagner  and  Manne  by  demonstrating  that  certain 
extra  constraints  involving  sums  and/or  differences  of  variables  from 
several  originsand/or  destination  constraints  are  in  fact  equivalent  in 
tre  transportation  setting  to  a  bound  on  a  parti *1  sum  of  variables  from 
a  single  origin  or  destination  constraint.  The  coefficient  stnjcture  of 
constraints  in  this  class  seems  upon  inspection  to  bear  no  relationship  to 
the  stapler  equivalent  forms,  ie  present  a  general  characterization  for 
icentifying  .  nese  constraints  and  give  a  procedure  for  generating  equivalent 
constraints  that  can  be  operated  on  by  Wagner's  or  Manne's  methods.  Ex¬ 
amples  of  orcbler  classes  in  which  the  more  general  constraint  structure 
appears  are  given  ir.  Section  4. 


2.  TRANSFORMABLE  CONSTRAINTS  AND  THEIR  EQUIVALENTS 


Tne  standard  transportation  problem  can  be  stated  as  follows: 


Minimize  Z  Z  ciixi 

M  jeN  J  J 

subject  to 

xij  1  0 

wnere  M  =  {l,...pn}  and  n  =  {l,...,n}  and  z  a^  =  z  b^ 

icM  jeN  J 

This  problem  is  frequently  presented  in  a  tableau  format  as  shown  in  Table 
1  with  a  sinqle  row  representing  each  origin  constraint  and  a  single  column 
representing  eacn  destination  constraint. 


r 


Our  "iain  assertion  is  that  any  singularly  constrained  transportation 

problen  whose  extra  constraint  is  of  the  form 

P:  :  :  -:ixij  <_  D 

i--M*  j  •-  0  -j  J 

where  M'  =  -  •  p;  for  some  p-:M  and  for  each  icM',  Ji  =  S  or  Ji  =  N-S 

.  .  J  +1  Ji  *  S 

ancl  -j-^-lifJ-j*  N-S  where  SCN 

can  Dfc  transferred  into  an  enlarged  transportation  problem.  To  give  an 
example  of  one  such  constraint,  consider  a  transportation  problem  with 
four  origins  and  five  destinations  as  depicted  in  Table  1.  The  associated 
supply  and  demand  values  are  indicated  on  the  right  and  lower  borders. 

A  specific  constraint  of  the  type  P  is 

X21  +  x23  *  x25  +  X31  +  x33  +  x35  'X42  'x44  -  9 
whose  non-zero  coefficients  appear  in  the  table  in  the  cells  corresponding 

to  the  appropriate  variables.  In  this  example 

S  =  ”-1,3,5;,  N-S  *  -.2,4;,  p  =  1 ,  Jg  =  {1,3,5;,  ={1,3,5;,  and  =  2,4; 


Demand 


C1 

32 

33 

°4 

D5 

°1 . 

°2 

+  1 

+1 

+  1 

°3 

+  1 

+1 

+  1 

°4 

-1 

-1 

7 

8 

2 

7 

6 

Supply 

8 

7 

10 

5 


Table  1 
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The  procedure  we  propose  for  transforming  such  a  constraint  into 
a  partial  sum  of  variables  in  a  single  rov  of  the  transportation  tableau 
creates  row  multiples  and  column  multiples  Kj  which  generate  a  linear 
combination  of  the  original  transportation  constraints.  Thereupon,  adding 
tnis  linear  combination  to  the  type  P  constraint  produces  a  constraint 
which  can  be  accomodated  by  the  methods  of  [5  ,10].  The  procedure  con¬ 
sists  of  a  single  iteration  of  the  following  instructions. 

Step  1.  Set  R1-  =  0  for  all  icM  and  Kj  =  0  for  all  jeN. 

Define  I ^  =  < i eM * :  =  S; 

]N-S  =  Ji  = 

Step  2.  If  Is  =  0  go  to  step  3.  Otherwise  set  R^  =  -1  for  all  i e  1^ 

Step  3.  If  N-S  =  t  go  to  step  4.  Otherwise  set  K.  =  1  for  all  jeN-S. 

J 

Step  4.  Multiply  the  ith  origin  and  destination  constraints  by 
their  corresponding  R ^  and  Kj  vali.es  and  sum  to  form  a 
linear  combination  of  the  standard  origin  and  destination 
constraints.  The  addition  of  the  resulting  equation  to 
the  constraint  of  type  P  produces  an  inequality  bounding  a 
sum  of  variables  in  row  p. 

To  illustrate,  we  apply  this  procedure  to  the  exam^e  in  Table  1, 

Step  1  sets  the  multiples  Rj ,  R2,  R3,  R4,  and  K-j ,  K2,  K^,  K^,  and  K& 
equal  to  zero  and  identifies  the  sets  I5  ={2,3}and  IN_^  =  {4}.  In  step  2, 

R2  and  R3  are  set  equal  to  -1.  In  step  3,  K2  and  K^  are  set  equal  to  +1. 
The  constraints  of  the  transportation  problem  associated  with  non-zero  R^ 
and  Kj  multiples  and  the  extra  constraint  are  shown  below  in  detached 
coefficient  form  with  the  non-zero  R.  and  Ki  multiples  indicated  on  the 
left. 
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:llx12x13x14x15x21x22x23x24x25x31x32x33x34x35x41x42x43x44x45  RrtS 


R2=(-l} 

V'-1* 
<2=(+1 ) 
•<4=(+U 

extra 


11111 


11111 


1  1 


1  1  -1 


Forming  the  linear  combination  described  in  step  4  and  adding  this  to  the  extra 
constraint  yields  +  x^  i  7,  as  shown  below  the  line.  Using  the 

techniques  described  by  Wagner  [10  ]  and  Manne  [  5,  p.  332],  tne  original 
problem  can  now  be  reformulated  as  an  equivalent  uncapacitated  transporta¬ 
tion  problem  with  one  additional  origin  and  one  additional  destination. 

The  validity  of  the  transformation  procedure  is  established  in  the 
following  theorem. 

Theorem:  The  addition  of  a  constraint  of  type  P  to  the  linear  combination 
of  the  standard  origin  and  destination  constraints  specified  by 
the  procedure  yields  an  equivalent  constraint  bounding  a  sum  of 
variables  in  row  p. 


Proof:  The  definitions  of  M* ,  1^,  IN_S»  and  give 


^  1  :  :  - ixi j =  $  1  xij  "  ‘  xij 

i cM*  j.;Ji  i cig  jeS  lcl;j_s  jeN-S 


Adding  tne  appropriate  linear  combination  of  the  transportation 
constraints  to  the  expression  on  the  right  above  yields 
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i.lls 

z  x 
jeS 

ij 

Z  z  Z  +  Z  l  X^j 

ie^-sjcN-S  J  iel$  jeN  J  ieM  jeN-S  J 

ilh 

jeN-S 

Xij 

-Z  l  x^  •  E  Z  X^ 

i cl,,  c  jeN-S  J  ieM  jeN-S  J 

if  J 

icM‘ 

Z 

jeN-S 

x-u 

+  u  L  xij 

ieM  jeN-S  J 

z 

u-:p> 

jeN-S 

xij 

=  1  Xpj  -  °' 

jeN-S  w 

where  D'  =  D  +  z  ♦  z  =  D  -  I  ♦  z  b.-. 

ieM  jeN  i  e.  I  ^  jcN-S 

This  completes  the  proof. 

3.  SOME  EXTENSIONS 

The  general  form  of  the  constraint  P  is  stated  as  a  "less  than  or 

equal"  constraint  type,  but  the  procedure  is  clearly  valid  for  both 

equality  and  "greater  than  or  equal"  constraints.  In  the  latter  case, 

an  inequality  of  the  form  z  xD^  >_  D'  will  be  produced,  which  may 

jeN-S 

be  made  suitable  for  transformation  by  Wagner's  and  Manne's  methods  by 

noting  that  it  is  equivalent  to  Z  xpj  <  ap  -  D'.  Additionally,  since 

jeS 

the  index  sets  M  and  N  are  arbitrarily  associated  with  origins  and  des¬ 
tinations  (or  rows  and  columns),  P  could  also  describe  a  class  of  con¬ 
straints  whose  members  are  equivalent  to  a  partial  sum  of  variables  in 
a  single  column. 

It  should  be  clear  that  the  5^  values  in  constraint  P  could  be  any 
value  so  long  as  all  ^  are  equal  in  absolute  value.  In  this  case  step  2 
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wojld  need  to  be  altered  to  set  the  appropriate  to  -k  and  in  step  3  the 
Kj  multiples  should  be  set  to  +k  where  jd^J=  k  for  icM 1 .  After  finding 
the  equivalent  constraint  in  step  4  both  sides  of  the  constraint  would 
then  be  multiplied  through  by  1/k  to  obtain  unity  coefficients  for  the 
variables.  If  the  original  right-hand  side  value  of  the  extra  constraint 
is  an  integer  multiple  of  j  ^Jand  if  the  a^  and  bj  in  the  original  trans¬ 
portation  problem  are  integer  valued,  then  the  optimal  solution  to  the 
equivalent  enlarged  transportation  problem  will  be  integer  valued.  Thus 
our  procedure  gives  a  convenient  way  to  solve  certain  integer  programming 
problems  with  the  computationally  efficient  transportation  algorithm. 

If  the  original  problem  includes  several  constraints  of  type  P  then 
the  procedure  could  be  applied  to  each  extra  constraint  independently. 

In  this  case  if  the  transformed  constraints  involve  disjoint  sets  of  var¬ 
iables  or  if  the  sets  of  variables  are  nested  in  the  same  node  constraint 
then  an  equivalent  larger  transportation  problem  can  be  derived  via  [5,10  ]. 

The  presence  of  additional  constraints  necessarily  raises  questions 
of  feasibility,  and  finding  an  equivalent  restriction  may  shed  light  on 
the  feasibility  of  the  original  problem.  In  particular,  if  the  new  con¬ 
straint  formed  in  step  4  results  in  all  zero  coefficients  and  a  non-neg¬ 
ative  right-hand  side  then  the  extra  constraint  is  redundant.  Also, 
cnecking  the  sign  of  the  coefficients  and  the  sign  of  the  right-hand  side 
(in  relation  to  the  direction  of  the  inequality)  may  indicate  whether  the 
oroblem  has  any  feasible  solution  whatsoever. 

When  constraints  can  be  identified  as  being  of  type  P  then  the 
equivalent  restriction  can  be  found  directly  by  the  formula 

:  xni  i  0  -  kz  a,-  +  ki  b^,  where  k  =  UJ  and 

jtN-:  ie I s  jcN-S  J  1  1 
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the  procedure  of  finding  the  appropriate  linear  combination  need  not  be 
explicitly  applied.  This  is,  of  course,  a  direct  consequence  of  the 
proof  of  the  theorem. 

4.  APPLICATIONS 


Constraints  of  type  P  appear  in  many  problems.  For  example,  consider 
a  problem  where  warehouses  ship  a  perishable  product  to  several  markets 
and  because  several  of  the  warehouses  are  physically  distant  from  some  of 
the  markets  a  special  container  must  be  provided  for  each  item  shipped 
along  certain  routes.  If  the  supply  of  these  containers  is  limited  then 
the  problem  can  be  modeled  as  a  standard  transportation  problem  with  the 
additional  container  constraint.  An  example  is  illustrated  in  Table  4 
where  the  +1  values  denote  the  routes  that  require  containers.  The  in¬ 
dicated  constraint 


M3 

H4 

W1 

W2 

+1 

+  1 

W3 

+1 

+  1 

Table  4. 


is  transformable  and  is  equivalent  to  a  constraint  on  the  shipments  from 
W1  to  M.j  and  M3. 

A  type  of  "balancing"  constraint  is  also  represented  by  the  character¬ 
ization  P.  Suppose  that  several  warehouses  can  supply  several  markets.  Given  de¬ 
mands  for  the  product  and  figures  for  warehouse  supply  a  transportation 
model  can  be  formulated  to  determine  the  minimum  cost  shipping  schedule. 

If  in  addition  management  requires  a  certain  degree  of  specialization,  two 
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constraints  may  be  introduced  to  balance  the  shipments  from  certain  ware¬ 
houses  to  certain  markets.  Tables  5a  and  5b  represent  two  different  bal¬ 
ancing  constraints  when  the  right-hand  side  of  both  constraints  is  zero. 
The  constraint  indicated  in  Table  5c  is  equivalent  to  that  in  5a,  and  the 
constraint  in  5d  is  equivalent  to  that  in  5b. (Since  the  constraints  in 
Table  5c  and  in  Table  5d  include  disjoint  sets  of  variables,  Wagner's 
procedure  can  be  applied  to  incorporate  both  constraints  in  a  larger 
uncapacitated  transportation  problem.) 


Warehouse 


Market 


Market 


1 

2 

3 

4 

1 

2 

3 

4 

A 

A 

+1 

+1 

B 

-1 

-1 

B 

-1 

-1 

C 

+1 

+1 

C 

Table  5a.  Table  5b. 


Warehouse 


Market 

12  3  4 


Market 

12  3  4 


Balancing  constraints  or  other  constraints  of  type  P  also  appear  in 
models  of  trading  between  nations  that  seek  to  maintain  favorable  balances 
of  trade  by  legal  restrictions  on  their  Imports  and  exports.  In  such  a 
model  the  total  flow  of  goods  from  one  nation  may  be  required  to  differ 
by  no  more  than  a  specified  amount  from  the  flow  of  goods  in  the  reverse 


direction. 


FOOTNOTES 


It  has  been  brought  to  the  authors'  attention  that  these  procedures 
were  developed  in  linear  programing  courses  by  A.  Charnes  and  W.W.  Cooper 
at  Carnegi e-Mellon,  Purdue,  Northwestern,  and  the  University  of  Texas. 
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